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ABSTRACT 

The thesis develops a theoretical model for predicting 
ship interactions that should be useful in a ship trajec- 
tory analysis. The method assumes that when ships are 
close together, the potential field is the main source of 
the interaction forces. The flows around the ships are 
represented by using half bodies of revolution generated 
with axial singularity distributions. A second axial 
distribution of doublets were included to correct for 
cross flows and thus maintain the rigid body boimdary 
condition. Following the determination of the fluid flows, 
the forces were calculated using Lagally’s theorem as ex- 
panded by W. E. Cummins. 

The principal result was a computerised model to pre- 
dict forces resulting from steady, parallel motion of 
bodies of revolution moving in an infinite ideal fluid. 

The program is capable of handling ship geometries in rec- 
tilinear motion. It was found that the theoretical model 
produced similar results when compared with model tests 
conducted by D, W. Taylor, however, those tests were re- 
ported as being crude. This emphasizes the need to obtain 
a set of more valid empirical data. 

Thesis Supervisors Martin A. Abkowitz 

Title: Professor of Ocean Engineering 
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CHAPTER I 



INTRODUCTION 



1 . 1 General 

Interactions between ships operating in proximity can 
be thought of as a modern day phenomena. There are no knov/n 
references to it in the literature until after the emer- 
gence of steam driven ships. Along with the increase in 
shipping in shallow and restricted waters, this innovation 
contributed to the higher frequency of collisions. These 
collisions provided the operator with his initial encounters 
with ship interaction forces. It was many years, however, 
before the courts recognized the existence of this strange 
force which contributed to many of the collisions. The 
acceptance of these forces came only after many eyev/itness 
accounts described strange ship behaviors and it was con- 
cluded that there must have been a force that was not 
present during normal operations. 

With the increased density of shipping in restricted 
waters, especially for channels, it is not surprising that 
collisions involving interactions were first noticed for 
the shallow water occurrences. The effect of shallow water 
has persisted as being important in describing the inter- 
action phenomena. Shallow water does influence the flow 
around the ship, however, for depths greater than twice the 

draft of the ship, the effect of the ocean bottom is 
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expected to be small. This is based upon the interaction 
between ships in deep water, forces have been experimentally 
determined as the ship separation was increased. 



With the increased speeds and maneuverability intro- 
duced by steam propulsion, interactions other than with 
ships were first noticed and recorded. These involved the 
forces experienced by the position of the ship in the 
channel. From these forces, rules of thumb were formed such 
as the ship "smells the bank". This particular rule indi- 
cates the mariner realized that the ship interacts with the 
channel bank so as to sheer away from it. It is in this 
area of ship maneuverability in restricted water that a 
great many empirical studies have been made. With the re- 
sults of Bernoulli, the mechanism could be explained. Even 
today the operators aware of the Venturi effect tend to 
associated all interactions and explain them consistent 
with this law. This may give a rough guide, but does not 
contribute to the understanding of the process. It may at 
times even lead an investigator to the v/rong conclusions. 

The forces evolving from ship interactions has long 
been a controversial subject. The main interest has stemmed 
from the cases which result from collisions. The literature 
is sprinkled with accounts of major collisions. Included 
are discussions where important figures of the day give 
conflicting explanations of the ships' behavior. The first 
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scientific investigations involving model tests occurred 
in 1909 when D. W. Taylor presented his findings to the 
Society of Naval Architects and Marine Engineers. The dis- 
cussions of the paper demonstrated the debate at the time 
about the causes of the forces. 

The study of interaction forces is important for at 
least four reasons. First, they contribute to collisions. 
The understanding of these forces is necessary to help pre- 
vent future occurrences. Second is for the control of 
ships moving in proximity to each other, as done by several 
Navys to accomplish refueling of ships at sea. The third 
implication is in the design of future ships and their 
operating procedures. The adverse effects of ship inter- 
actions may be reduced by modifying design parameters and 
effective design and use of controls. Finally, for those 
ships in operation, the present rules for navigation do 
not recognize the interaction effects. Future study may 
show the need to modify these rules so as to provide safer 
movement of the merchant fleet. 

With this perspective, the research effort here at- 
tempts to determine the extent of the past efforts to de- 
velop a generalized theory for ship interactions which 
v/ould be practical to the engineer. The view is toward a 
general model, in particular one which could be used as 
part of a generalized program for predicting ship trajecto- 
ries. 
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1 , 2 Ship Interaction 



1,2,1 General 

What is meant by interaction forces? Any force that is 
only observed v/hen one or more ships (or other objects such 
as wall, etc, ) are present can be classified as an inter- 
action force. The motion of a ship through the water mani- 
fests itself by waves and the wake formed and by a moving 
pressure disturbance in the water. The moving pressure 
field will be identified as the potential disturbance where 
the waves will be referred to as the free surface disturb- 
ance. 

The interactions between ships operating on the surface 
of the ocean result from the influence of both the free sur- 
face and potential disturbances, whereas each appears to 
create forces by different mechanisms and differ in their 
spatial distribution. The potential field generates a 
distribution of pressures in the fluid. The waves also 
generate a pressure field. However, this rapidly attenuates 
with depth. The waves also change the water surface and 
this introduces forces. If the frequency is low, the long 
duration of the force will cause significant ship motions. 

An example is in the case of following seas where direc- 
tional instability may occur and result in a large response 
of the ship. 

Very close to the ship, the potential field dominates 
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the fluid flow pattern, where much further away the waves 
produce the flow disturbances. The induced flows act on 
other bodies (ships) placed in the perturbed stream to 
create interaction forces. 



1.2.1 Taylor Results 

One of the first scientific experiments performed and 
reported on the interaction between ships was conducted by 
David V/. Taylor at the Model Basin in V/ashington, D, C, 

The results are reported in the 1909 Transaction of 
Society of Naval Architects and Marine Engineers, A good 
summary of the work can be found in Reference 2, Here 
attention is focused on the quantitative measurements. 

The results are presented for the two model trials. The 
first is with relatively thin models and the second is v/ith 
fuller models. For convenience, the model characteristics 
are also presented. In each case, the forces are reported 
for only one of the models as indicated. 
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FIGURE I 

TAYLOR'S RESULTS 
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The original results were non-dimensional ized by the 
resistance of the model. Taylor comments that the results 
may have inaccuracies due to problems encountered in towing 
the models on a parallel course. He also stated that the 
results for the thinner models were less consistant. It is 
difficult to say if the data was erratic or if the results 
were counter intuitive, thus indicating to Taylor that they 
might be erratic. V/hatever the meaning, both results are 
presented and evaluated using the proposed theoretical 
models. Appendix E contains calculations of the possible 
errors in the data that might have resulted from misalign- 
ment in the models. 

The original data was modified to develop the total 
forces and moments on the models. Also a change in the 
non-dimensionalization was made. To accomplish these 
changes, the force acting on the model was determined by 
multiplying the data by the resistance of the model, Taylor 
states that the models were towed at speed between two and 
three knots. A value of 2,5 knots was assumed. The 
total force on the model is determined by adding the forces 
that were measured at each end of the model. The moments 
were calculated about the midship position. The following 
curves depict Taylor's data after the indicated changes. 

The curves are drawn for the ships in various rela- 
tive positions. The abscissa is scaled in terms of the 
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FIGURE 2 



TAYLOR MODEL TESTS 
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model length. The attraction and repulsion forces are self- 
explanatory and are applicable for the model passing on 
the starboard and port side, The moments are depicted 
as being either bow in or bow out. This indicates the di- 
rection in which the moment is acting on the model. The 
Thin Ship data represents Model Number 834 as it passes 
Model Number 838 . The Full Ship data represents Model 
Number 858 passing Model Number 866. 

1,2,3 Newton Results 

Experiments were carried out at the Admiralty Experi- 
mental Works (ASW), Haslar, England, between 1946 and 1948 
to test out the feasibility of the replenishment of war- 
ships with fuel tankers. These results were reported at 
the first Symposium on Ship Maneuverability at the David 
Taylor Model Basin in Washington, D. C, in i 960 by R. N. 
Newton, The tests conducted included an experiment with 
constrained models of the Battleship H.M.S. KING GEORGE 
and the R.F.A. OLNA. (Ship A and B Respectively) 

A.E.W, used freely propelled and controlable models to 
study the behavior during maneuvering alongside and breaking 
away. The work also contains some data resulting from 
tests conducted at sea with ships. The results are presen- 
ted here for the constrained models. The force moment act- 
ing on each ship has been converted to a position amidship. 
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FIGURE 3 



AE.W. MODEL TESTS 
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1 , 3 Previous Analytic Models 



1,3.1 Ellipsoid 

Sir H. T. Havelock proposed a model of ship interac- 
tions based upon potential flow. As an example, he cited 
the ellipsoid. In the discussion to Reference 8, the solu- 
tions were given for two unequally sized ellipsoids moving 
parallel in the fluid. The forces were calculated using 
Lagally's theorem. 

The follov/ing singularity distribution describes an 
ellipsoid in a uniform flow. 

SM(X)=2*C*Vhx 

WHERE 

1/C = 4 e/(l-e2) - 2 log [ ^)/ (1-e) ] 

V = Speed in FT/SEC 
X = Distance from midbody position 
e= Ellipsoid eccentricity 

Appendix D contains sample distributions for the 
ellipsoid model. Also, the results for the sample case 
cited by Havelock are recalculated in Appendix F, 
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1.3.2 Rankine Ovoid 



The same model proposed by Havelock for use with the 
ellipsoid could be applied to any other body shape for 
which the singularity distributions are known. The Rankine 
Ovoid is one of the simplest bodies to generate from sin- 
gularities in the fluid. It is formed by placing two equal 
but opposite singularities in a uniform flow so that their 
axis is in line with the incoming flow and the source is 
upstream of the sink. Given the singularities, it is easy 
to determine the size of the body that is represented in 
the fluid. Silverstein in Reference l4 developed the force 
on a Rankine Ovoid. Appendix F contains the recomputation 
of a case demonstrated in his paper. Even though the dis- 
tribution is relatively simple and the body shape can be 
determined, it takes some effort to solve the inverse prob- 
lem. The desired singularity strengths and positions for 
a given shape requires the simultaneous solution of two 
non-linear equations. 

The following equations describe the Rankine Ovoid s 
4-, (Strength) (Dist) ( 1 / 2 ) 



V = 



(V2)2 = 



[ (L/2)2 - (Dist)^ ] ^ 

4. (Strength) (Dist) 



V J{B/2)^ + (Dist)^ ' 



\7HERE 

V = Speed in ft/sec. 

Dist = Distance from the modbody that singulari- 
ties are located. 



Strength = Singularity Strength 
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The solution of the simultaneous equations was accom- 
plished on the computer .using s system provided program 
called Zeroin. A copy of the instruction sheet for this 
program is included in Appendix H. A tabulation of Ran- 
kine Ovoid distributions is provided in Appendix D. 

The experience with the system routine that solved the 
simultaneous equations should help future investigations. 
The Zeroin Program is one of two that were available. This 
particular one was selected mainly because it offered an 
automatic feature which permitted retrieval of intermediate 
steps in the solution so that the program results could be 
verified. It was found that an approximate solution to 
the problem using the semi -infinite body would provide a 
good starting solution for the program. The routine would 
not converge to the problem solution, however, until the 
magnitude of the initial solution vector was reduced in 
size by two orders of magnitude. This prevented the prog- 
ram from jumping past the solution and attempting to find 
another solution at very large values of the variables. 

1.3.3 Slender Body Approximation 

Newman has derived the forces and moments acting on a 
slender body of revolution moving parallel to a wall. In 
his derivation, it was assumed that the fluid was ideal and 
of infinite extent, thus ignoring any free surface affects. 
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In order that the boundary condition at the wall be satis- 
file, an image singularity distribution v/as provided. This 
is analogous to two bodies of revolution moving parallel 
and abreast through an infinite ideal fluid. 

The approach used was to determine the singularity 
distribution that v/ould develop fluid boundaries representing 
the rigid body moving through the mediiim. The singularities 
were sized using a two dimensional model and then trans- 
formed to three dimensional by assuming that the body was 
slender (rQ<< L). The location of the singularities was 
determined by satisfying the linearized boundary condition 
on the body surface. It was pointed out that for the body 
close to the wall, the singularity distribution could be 
approximated by a general bipolar distribution. This re- 
sulted in the singularities being placed along the length 
of the body in such a v/ay as to be displaced from the center 
axis. This formulation requires the wall be close 
(Z/L « 1). The displacement of the singularities is des- 
cribed by the following formulas, where the displacement 
from the centerline axis is equal to (Zq - a). 




Once the singularity distributions are known, the 
forces acting on the body are calculated using Lagally's 
theory; It is easy to verify that the singularity distri- 
bution remains within the body and thus is amenable to a 
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Lagally force calculation. The following diagrams show the 
distributions for a sample ship form. 

The procedure developed the singularity distribution 
using a slender body assumption. The resulting force in- 
tegrals were then modified. An approximation to the inte- 
gral is based upon the assumption that the part of the in- 
tegrand v/hich represents the singularity as a function of 
length can be treated as a constant value. 

Nevmian's v;ork produced the following results for the 
probl em. 



V/HERE = Separation of the body axis from the wall 



Newman also presented a sample case to demonstrate the 
results. The body was assumed to be an ellipsoid des- 
cribed in the following form: 




= Radial dimension at position X along the 
body 



r. (x) = r., * (1-if 2) 



^VHERE 



r^ = Maximum radius of the body 
X = Longitudinal position measured from the 



midship station. 
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Substituting into the Force Equation, 

Pj = - %'rrfV‘r.^^L'' [QcU^-Qa!)] 

WHERE 

Qj^ is the Legender Function of the Second Kind. 
AND ^ = ^*4, 



AS 



*e MAX 



S-* i 

- -Z'tl-flf'fim L 



t 



The force integrals were evaluated for this case and 
the results presented as follows, 



FIGURE 4 
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In the concluding comments of the paper, the author 
states, ",.,.it seems likely that for bodies with full or 
blunt ends this [moment] integral will be dominated by the 
large values of S’(x) at the two ends, and thus the pitch 
moment will be positive (bow down) if the bow is blunt 
compared to the stern, and vice versa". In this comment, 
the pitch is analogous to the yaw in the previous defined 
orientations. He is predicting a ’bov/ out’ condition if 
the sterns are blunt compared to the bov/s. The author 
also notes that the force is always attractive and for 

geometrically affine bodies, having the same cross sec- 
tional distributions but different lengths, the force will 
be inversely proportional to the length whereas the moment 
will be independent of the length". From this result, 
the nondimens ional value of maximum force is inversely 
proportional to the cube of the L/B ratio. 
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FIGURE 5 
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FIGURE 6 
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CHAPTER II. 
FORCE MODEL 



2.1 Model Description 

The proposed model for predicting ship interaction 
forces uses the hypothesis proposed by Havelock that the 
potential field could explain the ship interaction pheno- 
mena. This approach seems reasonable for several reasons. 
First, the ship moving at the water's surface produces a 
flow field similar to that produced in an infinite fluid. 
Second, the potential field dominates the flow near the 
body. Pressures in this region fluctuate spatially and the 
changes find relief on the surface. The movement of this 
pressure field causes the wave system. Third, for most 
cases of interest, the ships are operating below their de- 
sign speed. This is particularly true for naval ships re- 
plenishing at sea. The fine entry of these ships causes 
minimal wave disturbances when at replenishment speeds. 

The potential flows without a free surface can be 
modeled using bodies moving in an infinite fluid with image 
systems to maintain the flow conditions on the waterline 
plane. This condition requires that there be no normal 
component of flow at the interface representing the undis- 
turbed v/ater level, A body of revolution moving horizon- 
tally in the fluid has a symmetry of flow and is used to 

model the ships. The maximum diameter of this body is taken 
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as the ship's beam because of the major influence this di- 
mension has. 

The force acting on the total body of revolution would 
represent twice the force experienced by the ship. The 
force results from the integration of pressures on (the 
lower) half of the body of revolution. The symmetry of 
flov/ will insure the same pressures and thus forces on the 
other half body. 

The major reason for selecting a body of revolution v/as 
the ability of using an axial distribution of singularities 
to describe the flov/, Hess et. al. have shown the compu- 
tations required for representing an arbitrary body in an 
ideal fluid by a surface distribution of singularities. 

The computational effort for each ship would be much greater 
for surface distribution and as such is not desirable for 
use in a ship trajectory model. 

In total, four models were developed and programmed. 

The first two incorporated the ellipsoid and Rankine Ovoid 
as discussed earlier. The third and final one involved a 
modified ellipsoid and slender body respectively, each with 
cross flow corrections. The first two used the derivations 
of Havelock and Silverstien. The third is only a special 
case of the fourth model. The slender body model will be 
discussed in more detail. Appendix I contains a logic flow 
chart of the computer program for this model. 
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The slender body model uses bodies of revolution to 
model the ships as described previously, the axial distri- 
bution of sources and sinks is obtained by assuming the 
body of revolution is slender. In two dimensions, the sin- 
gularities are sized using the following equations. The 
derivation for three dimensions are presented in Appendix C. 
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This axial distribution of singularities determines the 
body in a uniform flow. When another body comes near, the 
flow field is changed, so that it is no longer uniform. 

As the bodies approach each other, their affect in altering 
the flov/ field causes a deformation in the fluid boundary 
which describes the rigid bodies. To handle this problem 
an additional distribution of dipole singularities are 
placed along the axis, Thes- are oriented in the plus or 
minus direction on the body axis and are sized to counter 
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the effects of local cross flow on the bodies. The change 
in the flow in the axial direction is handled by perturbing 
the original singularity distribution. 

The boundary condition on each body gives no flov/ 
across the boundary which describes the rigid body. In our 
case, this is a body of revolution. At each location along 
the body there is a radial distance at which the velocity 
(radial) must be zero. If the cross flows were uniform, 
the motion of an infinite cylinder could be represented 
exactly in the fluid by placing the appropriate doublet 
distribution on the axis. The flow pattern around a slender 
body of revolution can be approximated locally by an infi- 
nite cylinder because the change in diameter would be small 
for a change in the axial direction. For the distances 
between bodies of interest, it was estimated that the uni- 
form flow assumption for any one segment of the body would 
be nearly correct. The dipoles were sized with flov/ of the 
value experienced on the body axis as a measure of the 
uniform flow strength. For the actual flow conditions, the 
boundary location is not calculated. It is assumed that 
the deviation from the exact position would be small. The 
dipole distribution is sized by the following rule. See 
Appendix C for derivation. 

CX ) = ' V(Ar) » /UfBA (y) 
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If the perturbed flov/ is known, the changes in the 
singularity distribution can be readily determined. For 
the case of two bodies near each other, the change in the 
singularitjr distribution in one of the bodies is reflected 
in a slight change in the perturbed flow field experienced 
by the others. V/hat is desired is a simultaneous solution 
to the problem. This is achieved by performing several 
iterations of the flow. Taking first one, than the other 
body, change its singularity distribution to counteract 
the perturbed flow, and then determine the resulting 
affect on the flov/ field at the other body. Changes in 
the perturbed flow were insignificant after about three to 
five iterations except when the bodies were etremely close 
and then it took only about ten iterations to achieve con- 
vergence to the solution. 

2.2 Interaction Force Calculations 

The interaction forces acting on a rigid body moving 
in an ideal fluid v/ithout viscosity can be determined if 
the fluid velocity field can be represented by a potential 
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function. In principle, the method for determining a force 
on any body would involve the calculation of the pressure 
on the body surface followed by the integration to obtain 
the total force. This was first accomplished by either 
Munk or Lagally for bodies generated with singularities. 
Lagally noted that for a body placed in a steady stream 
this integration took a form that would be analogous to a 
derivation which treats a force as acting between the sin- 
gularities. The apparent force between a pair of singulari- 
ties is determined in three dimensions by the follov/ing 
formula. 



The total force would be determined by adding all 
pair contributions. The force acting on a body which is 
represented by a continuous source distribution along its 
axis has the follov/ing integral form. 



The results of Lagally were expanded by Cummins to 
include the forces between doublet singularities for the 
general time dependent case. His derivation involved the 



F,2. "" 




Where a positive force indicates 
an attraction. 




Where [a.b] and [a',b'] indicate 
the line segments with singularities. 
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following assumptions. 



1, The velocity field is irrotational 
and has a velocity potential ? (x,y,z,t) 

2, If the body were not present, the 
stream would have a velocity potential fi, 
which v/e call the potential of the "undis- 
turbed stream" , 

3, There are no singularities of the 
undisturbed stream in the region occupied 
by the body. 

The boundary condition at the sur- 
face of the body is satisfied by super- 
imposing a system of singularities upon 
the undisturbed stream, such singularities 
falling v/ithin the region which the body 
would occupy. The potential of the system 
of singularities is designated by Then 



J = ^ + 

Cummins found that the forces and moments calculated 
were of the same form as expressed in the Lagally deriva- 
tion, but with additional terms. He expressed the Lagally 
force term as a product between local singularities and 
the total flow perturbation rather than as the sum of in- 
dividual contributions from each singularity. The two 
representations for the force are identical, Zucker gives 
a summary of the results and indicates the contribution 
by Landweber and Yih in substantiating Cummins' results 
and evaluating the moment term resulting from the change 
in the fluid flow with time. 

The force on the body will be as follows for Cummins' 



assumptions. 
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P = ^ V + F, Fz: ^ Fj 

' St 

Fi = f [Fj V + M, + K\3 

' dx 

WHERE 



"V" - Volume of the body 

Vo - Is the translational velocity of the body 

n = Is the location of the centroid of the 
' body relative to the origin of the moving 
system, 

F' The force if the body were not rotating 
and the undisturbed stream were steady 
(Lagally Force) 

-Force due to change of the flow with time 



Fj-Force arising from body rotations 



l^i -Moment if the dlow were steady (Lagally) 



H ~ Moment due to the flow changing with 
time 

' Moment effect due to rotation of the 
body 

FOR SITUATIONS DESCRIBED BY ONLY SOURCES AND DOUBLETS, 
F^ AND ARE ZERO. 
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FOR A SOURCE OP STRENGTH 



WHERE defines the fluid velocity in the 

undisturbed stream 

n - 

^ - y F, - 

FOR A DOUBLET QZ MOMENT 

-fVf (/'V) ^ 

- ^-TT f (/I 

Si-' V^x 

NOT EVAIiUATED BY CUMMINS ARE THE MOMENT DUE TO THE 
FLOW CHANGING WITH TIME. 



' H /^(+f) i 

wHe/tc: n= /-yiA//*<<t> tx^fcrft v/tcro< 



This term was evaluated by Landwebber and Yih 
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The force calculations programraed on the computer 



used at least one of the forms discussed previously. As 
more complicated and detailed models were developed the 
form of the force calculation was selected to improve com- 
putational efficiency. For the Rankine Ovoid and ellipsoid 
which do not contain cross flow corrections, the simplest 
statement of Lagally's force was used. V/hen it became 
necessary to introduce dipoles to correct for the cross 
flow, the final form of Cummins was used. The final 
slender body model takes into account the terms that are 
represented by Fj^ and Mj. . This represents the non-zero 
terms of the interaction force for steady parallel motion 
of the bodies. 

The above equations relate to a three dimensional 
flow problem. The two dimensional case is similar with 
the Lagally force being of the following form: 



Havelock developed in Reference 7 the force and moments 
acting on a two dimensional body by using contour Integra 
tion. 

FOR SOURCE DISTRIBUTION 



F;. r 2D-P Vr;, 



X ' , -2, C,x,-x.) 
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WHERE 



1. Suffix "S" refers to the given distribution 

in the liquid, and "r" to the image system 
within the surface of the body. 

2. Summation extends over the external and inter- 

nal sources taken in pairs. 



FOR DIPOLE DISTRIBUTION 



X - - ■‘It- f 2 Mj ) / , 

(h,s) ' Ky5 



Y' ^ " -'mS '3SrJ A 3 

(r.j) fr.O / ^*^5 



Where: — the doublet is of moment and makes 

an angle o{ with the x axis (ship). 

And; Ots being the angle between x' and the vector 
drawn from the doublet (r) to the doublet (s). 

The total moments results from moments of forces already 
determined plus an additional term. 




a-n-p 
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2. 3 Computation 



The earlier discussion of the slender body model gave 
the general steps involved for determining the interaction 
forces. In this procedure there were two major areas 
requiring considerable effort. First, it was necessary that 
derivatives of the flow potentials be derived for the force 
calculations. These potentials and their derivatives were 
defined in terms of a moving axis system of the ship in 
which the forces were acting. The necessary transformations 
and derivations are developed in Appendix B, 

The second area involved the determination of the total 
force acting on the ship. The slender body was defined by 
a continuous distribution of singularities. The resulting 
flows in the fluid v/ere determined by integrating the con- 
tributions from the singularity distribution. Also the 
forces were determined by integrating the individual force 
contributions. In both cases, the numerical integration 
used values of the integrand at the same location on the 
body axis. These coincided with the stations at which 
the sectional area curve had been defined. The integration 
was performed in each case using Simpson's Rule with an 
odd number of stations (21), evenly spaced along the body 
axis. 
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CHAPTER III. 



COMPARISON OF RESULTS 
3.1 Taylor Model Tests 

The Rankine Ovoid and ellipsoid models were developed 
and checked with the available results then used to model 
the forces and moments reported by Taylor. Figure 7 de- 
picts the results. These predictions differ in two v/ays 
with Taylor, First, they do not predict the secondary re- 
versal in the moment that is evident at longitudinal se- 
parations of ,6*h. Secondly, the largest moment is pre- 
dicted for longitudinal separations above 1.0 L. This 
variation can be attributed to the deformation of the boun- 
daries of the bodies as they approach each othar. If the 
boundary was maintained in the original position, the flow 
would have been restricted, thus increasing the fluid 
velocity and contributing to a lower pressure on the bodies. 
This reduction in pressure would contribute a moment that 
would cause a reversal in sign. The inability of the pre- 
vious theoretical models to explain the reversal in sign 
of the moment curve, led to the development of a model which 
incorporated cross flow corrections. The correction for 
cross flow indicates the reversal in the moment can be 
modeled. However, the magnitudes were not in full agree- 
ment with those obtained by Taylor. Computations were made 
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FIGURE 7 
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at higher speeds with no variation in the results. This 
was expected, the plot involved the non-dimensional moment 
coefficient. For closer distances, the predicted ampli- 
tude of the moments increased only slightly, (The model 
tests were conducted at a very small separation, ) 

The dimensions of the models used in Taylor's experi- 
ments were inputs to the slender body model. Calculations 
were made for the same cases as described by Taylor, The 
following curves show the results plotted for comparison 
purposes. 

Additional curves are drawn which show the forces pre- 
dicted for the Taylor full bodied models, as the lateral 
separation is increased. The reversal in moment diminishes 
as the ships separate because there is no longer a restric- 
tion in the flov/ to cause a pressure drop between the ships. 
Beyond a distance of 2 beams, the results agreed with those 
obtained v/ith a model which did not incorporate cross flow 
corrections. Only beyond a separation of 4 beam widths 
interaction effects were negligable. 

The model tests conducted by Taylor, and the slender 
body predictions differ in several ways. Some of the more 
significant assumptions used for the slender body weres 

1. The potential field was assumed most 
important (No free surface) 

2, The ships were represented by a body 
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FIGURE 8 
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FIGURE 9 
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of revolution. (Lessen effective beam) 

3. Squat and sinkage were not incorpo- 
rated , 

Other differences in the results may be attributed to either 
asemmetry in the models or inaccuracies in the measurement 
techniques. The apparatus used by Taylor was not indicated. 
V/hen the tests were conducted, the resistance measurements 
used a system of weights, pulleys and levers. Taylor even 
commented on the problem of conducting the model tests. 
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3. 2 Comparison of Results with Newman 



In order that a comparison be made between the model 
developed and the results of Newman. The sample problem 
in Newman’s paper was used. The sectional area curve for 
the ships moving in parallel were determined using the 
radius formula presented earlier. Two ships moving in 
parallel were defined with the same geometry. The initial 
calculation assumed a length of 600 feet, a beam of 100 
feet and a speed of 10 knots. The lateral displacement 
was determined and plotted according to the non-dimensiona- 
lization used by Newman. 

Note that Zq/tq is equivalent to the separation dis- 
tance between the axis of the two bodies divided by their 
beam dimension. Other beam to to length ratios were cal- 
culated by varying the length of the bodies while main- 
taining the beams constant. These included the L/B of 6, 
10, 20, 6o. The total force on the body of revolution 
determined using the model was non-dimensionalized by the 
maximum value determined by Newman. The following graphs 
depict the results developed in this thesis. 

The first graph shows the results for a L/B of 6 for 
the various theoretical models. Calculations were then 
made for the resulting forces with and without the cross- 
flow correction which used the dipole distribution. It 
was expected that the results of Nev/man and the proposed 
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FIGURE 10 
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model would approach each other in the limit of high L/B 
ratios and for very slender bodies because the singularity 
distributions would be nearly identical and lie very close 
to that used by Newman. 

The next graph shows the change in the non-dimensiona- 
lized maximum force as a function of the length to beam 
ratio. Note that the singularity model approaches the New- 
man result in the limit. On the graph, the error of the 
computation is indicated. The solution with the dipole 
correction seems to converge initially and then for very 
high ratios diverge. This tends to indicate that the di- 
pole correction may be useful for all ratios of interest. 

The error introduced during integration was easy to 
approximate for the axial distribution of singularities 
because the fourth derivative of the integrand should be 
exactly zero. However, for the case with the dipole cross 
flov/ correction this may not be the case and it is expected 
that the introduced error would be larger. This may ex- 
plain the divergence of the solution for high L/B (large 
lengths). See Appendix E. for derivation of error estim.ate. 

Another graph shows more directly the variation in 
the results for the various length to beam ratios. In this 
graph the shaded area indicates the effect of the dipole 
correction. 

The last graph is a compilation from various sources 

of interaction forces. They are plotted in the same 
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FIGURE 12 
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non-dimensional terms. For the ships, the forces were 
doubled to approximate the forces that would be observed 
between bodies of revolution. The non-dimensionalization 
with respect to the maximum force predicted by Newman is 
consistant with the other graphs. Plotted also are the 
solutions predicted for a body of revolution with L/B 
ratios of 6 and 10 with the geometry used by Newman. The 
plot indicates that the model tests agree more closely to 
the results predicted by the proposed model with the results 
of the tests by Newton showing the greatest variation. If 
the results had been mislabeled the results would have been 
more consistant v;ith the proposed model predictions. 

The basic differences in the proposed model and that 
of Newman are, first, the surface condition on the body is 
accounted for by displacing the singularities from the 
centerline in the case of Newman, v/here the proposed model 
uses a dipole distribution. It is obvious that the pro- 
posed model will satisfy the surface condition better v/hen 
the induced velocity is fairly uniform over the body pro- 
file. From the data presented, the relatively small L/B 
and the separation distances of interest justify the use 
of dipole corrections. A second difference in the models 
is in the evaluation of the integral solutions, Newman 
uses an approximation to the solution which emphasizes 
the need to restrict the solution to the slender bodies, 
whereas this model uses numerical integration. 
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CHAPTER IV 



CONCLUSIONS 



The major effort of the thesis has been to provide a 
theoretical model for the prediction of forces between 
ships. A computer program has been developed which carries 
out the computations for ships in rectilinear motion using 
the potential flow model described. The ellipsoid model 
proposed by Havelock was unable to predict the form of the 
results as indicated by empirical model tests. These 
tests demonstrated a reversal in the moment force as the 
ship separation decreased. This led to the development of 
the slender body model with cross flow corrections. The 
programs which were developed for ships in rectilinear mo- 
tion ignore time dependent forces. The data presented 
here represents the interactions of ships in deep water, 
on parallel courses and with each ship at the same speed. 
Orientation, other than parallel, would have resulted in 
a change in the relative positions of the ships with time 
and would involve other force terms as derived by Cummins, 
but were not programmed here. 

The investigation of past efforts indicate a lack of 
good model or full scale data to test a theoretical result. 
The differences experienced in the analysis of Taylor's 
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results could not be fully explained without more empirical 
data. This can only be obtained with ship models because 
of the parallel course requirements. 

In meeting the objectives of the thesis, the axial 
distribution of singularities and the force calculation 
techniques provide additional inputs to a general trajec- 
tory model. When the ship approaches other bodies or 
boundaries in the fluid, this work will help to predict 
the interactive forces so that more inclusive trajectory 
calculations can be made. This is applicable to the deter- 
mination of the factors involved in ship collisions, the 
prediction of the control parameters for at sea replenish- 
ments, or for future ship design characteristics. 

The performance of model tests may indicate results 
that differ from Taylor's. The design of model experi- 
ments must take into account the large forces introduced 
by any misalignment of the models. Also, the wall effects 
may introduce erroneous readings. The computer program 
has been written to include more than two ships in proxi- 
mity. It could be used to estimate the wall effects in 
the tow tank experiments by introducing the first few image 
distributions to approximate the wall boundary conditions. 
The following figure shows the orientation. It will be 
important that some tests be repeated with the models in 
a reverse geometry to insure accurate readings. This would 
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give a measure of the wall effects and asymmetries in the 
models. 

The program can he extended to cover shallow water 
effects. For very shallow water, a tv/o dimensional model 
could be developed by only a few changes in the present 
program. These would involve the change in sizing the 
singularity distributions and in the Lagally force calcu- 
lation. For other shallow water conditions, the boundary 
conditions could be satisfied by a series of image distri 
butions as shown in Figure IS. 
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FIGURE 15 
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CHAPTER V. 



RECOMVIENDATIONS 



The v/ork accomplished in the thesis points out the 
need for continued study of the interaction between ships. 
The following list indicates possible areas for future in- 
vestigation. 

1. Perform model experiments to obtain data to 
verify Taylor's results, including shallow water effects 
as described. 

2. Include remaining force terms -in order to calculate 
for unsteady ship motions. 

3. Incorporate program in trajectory model. 

4. Determine stable positions alongside, and test 
with operating experience. 

5. Simulate refueling maneuvers and compare with 
operator experience. 

6. As a final test, attempt to reconstruct published 
collisions involving ship interactions. 

7. Investigate feasibility of extending the theoreti- 
cal model to include free surface condition by distributing 
singularities as used by Havelock and developing the re- 
quired theorems to evaluate the interaction forces as accom- 
plished by R. D, Zucker for vortex sheets. 
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APPENDIX A 



SSAC (I,J) 

SAC (I,J) 
DSAC (I,J) 

DISPL(I) 
LENGT(I) 
BEAFi (I) 
VELOC (I) 
CP (I) 

CM (I) 

DIST (I) 

S (I,J) 

SM (I,J) 

SD (I,J) 

NN (I) 

KK (I) 

X (I) 

Y (I) 

Cl (I) 



= Nondimensional sectional area at Station 
J of Ship I 

= Sectional area at Station J of Ship I 

= Derivative of sectional area at Station J 
of Ship I 

= Displacement of Ship I (Not used) 

= Length of Ship I 

= Beam of Ship I 

= Velocity (retilinear) of Ship I 

= Prismatic coefficient of Ship I 

= Midship coefficent of Ship I 

= Distance from midship where R. 0, singu- 
larities are placed to generate Ship I 

= Initial (no interaction) source strength 
at position J on the axis of Ship I 

= Source strength at position J on the axis 
of Ship I 

= Dipole strength at position J on the axis 
of Ship I 

= Number of stations used to input non- 

dimensional sectional area curve and to 
perform calculations 

= Index for Ship I to indicate model used 
to calculate forces (not used) 

= X position of Ship I in the fluid 

= Y position of Ship I in the fluid 

= Heading angle of Ship I in the fluid 
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PI 

RON 

M 

MT/1 

N 

XFl(J) = 

YFl(J) = 

RFI(J) = 
XP2 
J 
L 

X1(J) = 

Y1(L) = 

I 

K 

APX(I) = 

APY(I) = 
AM(I) = 
ACM(I) = 



The constant 

Fluid density 

Index, (not used) 

Index (not used) 

Number of ships interacting 
(Program can handle up to 5 v/ith space allo- 
cation) 

Component of force in the X direction at 
Station J 

Component of force in the Y direction at 
Station J 

Moment at Station J on Ship I 

Component of force in the X direction of Ship 2 

Index for indicating position along Ship I 

Index for indicating position along Ship K 

Longitudinal distance to Station J on Ship I 

Longitudinal distance to Station L on Ship K 

Index to indicate first ship of the pair for 
force calculation 

Index to indicate second ship of the pair for 
force calculation 

Total force on Ship (I) in X direction (Out- 
put on non-dimensionalized ED-form) 

Total force on Ship (I) in Y direction 

Total moment on Ship I 

Total non-dimensionalized moment for Ship I 



62 



QR 


(I, J) 




Total flow at Position J of Ship I (not 
used) 


QA 


(I,J) 




Total axial flow at Position J of Ship I 


QN 


(I,J) 




Total transverse flow at Position J of 
Ship I 


DQA 


(I, J) 




Total axial contribution of flow gradient 
at location J of Ship I 


DQN 


(I, J) 




Total transverse contribution of flov/ 
gradient at Location J of Ship I 


QAZ 


(L) 


rr 


Axial flov/ contribution at J on Ship I 
due to singularity at L on Ship K 


QNZ 


(L) 


= 


Transverse flov/ contribution 


DGNZ 


(J) 


=r 


Transverse gradient contribution 


QRZ 


(J) 




Flow contribution 


DQAZ 


(J) 


=r 


Axial gradient contribution 


DXI 


(J) 


=r 


Station spacing on Ship I between Station 
J and J+1 


DXK 


(L) 




Station spacing on Ship K between L and 
L+1 


AEX 


(I) 


=r 


Focal distance of ellipse measured from 
midship position 
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APPENDIX B 



Derivation of Axis System and 
Transformations and Potential derivatives 
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MOVING AXIS SYSTEM 
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DEFINITION OF AXIS SYSTEMS 
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CONVENTION FOR DETERMING -RELATIVE 

DISTANCES. BETWEEN SINGULARITIES IN 
THE FORCE CALCULATIONS OF SHIP 1 
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RESOLUTION OF FORCE ON SHIP 1 
INTO AXIAL AND NORMAL COMPONENTS 
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APPENDIX D 



SINGUT.ARITY DISTRIBUTIONS 

The following graphs depict some of the representative 
singularity distributions. 

Model of an ellipse with L/B=6, L= 600 ft. , V= 10 kts. 
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The modi;f'ied ellipsoid has three singularity distri- 
when placed moving parallel to a wall with 
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The distributions used to model the Taylor full ship 
results are as follows. The situation when the models 
are abeam with a 3.9 ft. separation is represented (end 
corrections included). 
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APPENDIX F 



Presentation of Ellipsoid and 
Rankine Ovoid Results (Compare with original work) 
The results are the same except for 
a scale change. Havelock's scale should be divi- 
ded by 10, where Silverstien' s should be divided 
by (,4'rr')*^ or 158, 
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APPENDIX G 



COMPUTATION E?^ICIE^gY 

The programs written without cross flow corrections 
take only a fraction of a second to evaluate interaction 
forces for each individual ship position. The slender 
body model with cross flow corrections takes various times 
for execution depending on the number of iterations of 
the flows required. For the samples run, tv/o ships with 
twenty one stations each, took three seconds for five 
iterations increased to ten, the time doubled to six sec- 
onds for each ship position calculated. 
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Applications Progrcira Series 
AP-21 revision 2 



APPENDIX H 



I NAME : 

DESCRIPTION : 



METHOD: 



USAGE : 



MASSACHUSETTS INSTITUTE OF TECHNOLOGY 
INFORMATION PROCESSING CENTER 

September 28, 1973 

ZEROIN, a double precision FORTRAN IV subroutine 
to solve a system of simultaneous, nonlinear equa- 
tions 

ZEROIN computes a solution vector X of a set of N 
simultaneous, nonlinear equations: 

F^(X)=0, 



F^(X)=0. 

X is obtained by an iterative method; 

^(k+i) ^ ^(k)^ ^(k) 

where the vector *F (X , with 

(k) 

F(X ) representing the vector of function 

(k) (k) 

values at the point X and J representing 

(k) 

an approximation to the Jacobian at X . The 
following test is used to determine convergence: 

I 11*1 EPS* I I I * 

where EPS is a user-supplied tolerance. All real 
arithmetic in ZEROIN is double-precision. 

The call for ZEROIN is : 

CALL ZEROIN (N,X,R, EPS, ICV,EVAL,IRITE) 

Input values are; 

N the order of the system of equations, 

2 £ N < 20. 

X a one-dimensional array initially con- 

taining an estimate of the solution 
vector. 
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Page 2 



R 

EPS 
I RITE 



EVAL 

Output values 
X 
R 

EPS 

ICV 



an estimate for | j solution vector~X| | 

a tolerance to be used in the test 
for convergence, usually 0.5E-12. 

a specification for intermediate 
printing. If IRITE = 0, ZEROIII does 
no printing. If IRITE = 1, ZEROIN 
prints, at each iteration, the current 
value of R, the current approximation 
to the solution, the number of times 
EVAL (below) has been called, and the 
current function values F(X). * The 
printed output from the Tinal itera- 
tion contains, in addition, the output 
value of ICV. 

is the name of a user-supplied 
SUBROUTXUE subprogram that evaluates 
the functions F at each iteration. 

The call for EVAL is CALL EVAL(X,F) 
where X is the input to EVAL, and F 
contains the function 'valuers corres- 
ponding to X. X and F must be singly- 
dimensioned, double-precision vectors. 

will be: 

the solution vector. 

I I R /I II" 

where is the solution vector. 

When ZEROIN terminates normally, the out- 
put value of R v;ill be less than the in- 
itial value of EPS. In 'case of abnormal 
termination, R can be used to determine 
what degree of convergence was obtained 
in the output value of X. 

I I I ^ where F(X) is the vector of 

function values corresponding to the 
solution vector X. 



one of the following indications: 



ICV = 


0 


if convergence was obtained, 
1 |F(X) I P > 0.5E-10 


but 


ICV = 


1 


if convergence was obtained, 
1 l£(X) 1 P £ 0.5E-10 


and 


ICV » 


2 


if convergence was not obtained 
within M iterations, where 



M « MIN(10m,50) . 
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CODING 

CONSIDERATIONS : 



AVAILABILITY; 

I 

I CORE 

I REQUIREMENT: 
REFERENCES; 



AUTHOR: 



Appropriate declarations are: 

REALMS X(N) , R, EPS, EVAL 

INTEGER*4 N, ICV, IRITE 

EXTERNAL EVAL 

1. The input value of R may be chosen to be 

I |x| 1/10.0, where X is the input vector of es- 
timates, if the .solutions are all about the 
same magnitude. 

2. The output valuii of R Can be used to determine 
degree of convergence in the case ICV=2 . 

3. ZEROIN is more likely to converge if the equa- 
tions are scaled so all solutions are of the 
same magnitude, and if the initial estimates 
are good. Poor initial estimates may produce 
divide check interrupts during execution of 
ZEROIN. For details on the choice of initial 
estimates, refer to the first reference below. 

ZEROIN is available in load module form in the li- 
brary SYS5.MATHLIB.SUBR and in source fom in 
SYS 5. r^ATHLIB. SOURCE. For further information about 
the use of these libraries, refer toAP-60 or con- 
sult the Programming Assistants in Room 39-219. 

On the 370/165 ZEROIN uses 9,730 bytes of core stor- 
age (hexadecimal 2602) . 

Mancino, O.G. ** Resolution by Iteration of Some Non- 
linear Systems,” Journal of the Association for 
Computing Machinery , Vol. 14, No. 2 fApril, 1967) , 
pages 341-350. 

Robinson, S.M. ”Interpolative Solution of Systems 
of Nonlinear Equations," SI7U4 Journal on N umerical 
Analysis , Vol. 3, No. 4 (December, 1966), page 650. 

ZEROIN was written by G.W. Westley, Computing Tech- 
nology Center, Union Carbide Corporation, Nuclear 
Division, Oak Ridge, Tennessee. 
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Computer Plov/ Chart 
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RELATIONSHIP BETWEEN PROGRAMS 
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Computer Source Listing 
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